Two-component sandpile model : self-organized criticality of the second kind 
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Two-component sandpile models are investigated numerically and theoretically. Monte Cairo 
simulations are performed to show that probability distribution functions of avalanche size and 
lifetime obey power laws whose exponents are approximately equal to 1.5 and 2.0 and the system 
exhibits SOC. A mean-field theory is developed to discuss the essence of the processes. We find that 
two-component models approach a steady critical state belonging to a different universality class 
from that of one-component models. Conservation of two kinds of sands at local toppling causes 
an infinite number of stable states which substitute for artificial boundary dissipation. Among two 
control parameters appearing in one-component models, therefore, a rate constant of dissipation is 
removed in two-component models. It is concluded that the more conserved quantities result in the 
less control parameters and a novel class of SOC. 

PACS numbers: 05.40.-a, 45.70.Cc, 05.65+b 
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According to a vast amount of studies on fractals, 
1// noise, and so on, it has been reported that a 
lot of nonequiliblium systems in nature exhibit power 
laws. However it has not been deeply understood why a 
huge variety of distributions spontaneously go to power 
laws. Self-organized criticality(SOC) proposed by Bak et 
al. 0,0, SE1 has provided one reasonable understanding 
of the emergence of power laws through sandpile models. 
The models are cellular automata evolving in threshold 
dynamics. The domino effect induces a sequence of top- 
plings called an avalanche, whose distributions of mag- 
nitude and lifetime obey power laws without tuning any 
parameters. A number of works on SOC have been done 
numerically. Several analytical results have also been ob- 
tained. DharH has found some exact solutions in 
abelian sandpile models and has shed light on the stat- 
ics of the models, such as the number of total recurrent 
states and hight correlation functions. However, the dy- 
namics of sandpile models is still obscure, that is, the 
power-law distributions of avalanches still have been un- 
able to be derived analytically. 

Though sandpile models seem to have no controlpa- 
rameters by appearances, Vespegnani and Zapperi[J El 
have found that the models do have hidden control pa- 
rameters: a slow dissipation rate and a slower addition 
rate of sands. To explain this, they have introduced a 
mean-field theory of rate equations, 



toppling rule, it is concluded that the models become 
critical in the double limit 
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Pk = f K {Pa, Pc, Ps), K = a,C,S 



(1) 



where p a , p c , and p s are the probability densities of ac- 
tive, critical, and stable states, respectively. Imposing 
the conservation law of the number of sands in a local 



0, p a = S/e 0, 



(2) 
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where e is a dissipation rate and 6 an addition rate. In 
other words, p a is an order parameter and e, S are con- 
trol parameters. SOC is achieved by rough tuning of 
the " unapparent" parameters e, S around zero. Most of 
SOC modelsH [lfj belong to this "universality" class. 
We call this class SOC of the first kind. It seems inter- 
esting whether other kind of SOC exists, or equivalently, 
whether the condition of the double limit J2J are strictly 
the necessary condition of SOC? In this letter, we give 
an answer to these questions. 

As mentioned above, the conservation law of the lo- 
cal toppling rule plays a key role for SOC. Tsuchiya and 
Katori 11,] have proved rigorously that SOC breaks down 
in abelian sandpile models with nonconservative toppling 
of sands. Here, we pay attention to the number of con- 
servation laws. Does anything new happen if the number 
of conservation laws is increased? How about the robust- 
ness of SOC? To answer these questions, we consider two- 
component sandpile models which deal with two kinds of 
sands. This means the models have two kinds of con- 
servation laws. Firstly, numerical simulations are carried 
out to check power-law behavior. Secondly, a mean-field 
theory is developed to examine the essence of the pro- 
cesses. 

A two-component sandpile model is a cellular automa- 
ton defined on a regular lattice. A pair of two non- 
negative integers h(x) = is assigned to each site 
x on the lattice, where i,j{> 0) denote the numbers of 
sand A and B. At each time step, one unit of sand A or 
B is added to a randomly chosen site at a relative ratio 
(0 <)tab{< 1) of sands A to sands B. Several types of 
toppling rule can be adopted. For example, a toppling 
occurs when (a) i > ith or j > j t h, (b) i > ith o,nd 
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FIG. 1: Stable and unstable states in the rule (b). 



j > jth, or (c) i + j > k t h > 0) where i th , jth, and 
k t h are certain threshold values. If the number of sands 
at one site reaches threshold, sands A and B on the site 
topple one by one to randomly chosen nearest neighbor 
sites until the number of sands at the site becomes less 
than threshold. In these rules, a sand A and a sand B 
topple jointly. For example, if a site has n grains of sands 
A and no sands B, and then m (< n) new grains of sands 
B are added, it ends up with n — m grains of sands A 
and no sands B. These rules satisfy two kinds of conser- 
vation laws of local toppling. A series of topplings (an 
avalanche) continue unless the numbers of sands at all 
sites become less than threshold. Among three rules of 
toppling, the and rule (b) is the most interesting because 
of the presence of an infinite number of stable states. In 
the rule (b), stable states h s and unstable states h u are 
defined by 

h s = | < i < i th or < j < j th }, (3) 

h u = \i>hh and j > jth}. (4) 

and there are infinite stable and unstable states as il- 
lustrated in Fig. ^ These stable states can absorb all 
the added sands. This means avalanches can stop with- 
out introducing artificial boundary dissipation, which is 
indispensable in one-component models. Rules (a) and 
(c) are not capable of removing boundary dissipation be- 
cause the number of stable states is finite under these 
rules. Hereafter, we consider only the rule (b) and treat 
threshold ith = jth = 1 and periodic boundary condi- 
tions. 

Distribution functions of avalanche size S and lifetime 
T are calculated numerically, where S is defined by the 
total number of topplings in one avalanche and T by the 
total time steps of simultaneous updates of topplings. 
Cumulative distribution functions (CDFs) of avalanche 



size D(S) and lifetime D(T) in 1— dimensional lattices are 
shown in Figs. [21 We find that the CDFs obey power laws 
D(S) ~ S~ T and D(T) ~ T~ a where the exponents are 
approximately t ~ 0.5 and a ~ 1.0. Equivalently, prob- 
ability distribution functions (PDFs) are represented as 
P(S) ~ S~( 1+T ^> and P(T) ~ S~( 1+a l These exponents 
1 + r ~ 1.5 and 1 + a ~ 2.0 are very close to the criti- 
cal exponents of branching processes [l^- We carried out 
simulations in 2— and 3— dimension and obtained almost 
the same CDFs as those in 1— dimension. These results 
imply the model has meanfield-like characteristics^^. It 
is also found that varying tab from 0.5 to 0.1 hardly 
affects power-law tails. This means the model ends up 
with going to SOC states robustly in the long time limit 
without tuning tab- It becomes evident that the two- 
component sandpile models exhibits SOC. It should be 
noted that because of the absence of boundary dissipa- 
tion, the power-law tails become extended infinitely as 
the number of added sands N increases even in a finite 
system. 

Next, we construct a mean- field theory of two- 
component sandpile models and compare with that of 
one-component models^, @- We introduce the sets of 
probability densities of stable states X and unstable 
states Z as 

X = {X {id) | < i < 1 or < j < 1}, (5) 
Z = {%,) |i>landi>l}. (6) 

Rate equations for X and Z are given by an infinite 
system of first-order nonlinear ordinary differential equa- 
tions. However, it is difficult to handle the infinite de- 
grees of freedom. In order to avoid the difficulty, we 
define the following reduced six variables 

X* = {X , Xa, X B }, Z* = {Z 0l Za, Z B } (7) 

with 



Xq = X, 
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x A = ^x (ifih z A = ^2^2z (iJ) , 

i=l i=2 j=l 

oo oo oo 

X B = Y^ X(o,j) , z b = /ZY1 Z ^i) ■ 

3 = 1 i=l j=2 

For obtaining a closed set of equations for X* and Z* , 
we assume Poisson distributions for Xu t0 \ and -X7o,j) as 



X, 



(i,o) _ Ma 



i-l 



X, 



X A (< - 1) 

3-1 



(0,3) _ Ms 



Xb (j - 1) 



e~ MA (i > 1), 



(j > 1) 



(8) 
(9) 



where (1a = $At> Ms = $Bt are mean values of each 
Poisson distribution at time t and 5a, 6b rate constants 
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FIG. 2: Double logarithmic plots of CDFs of avalanche size (a) and lifetime (b). The lattice size L = 10 4 and the number of 
added sands N — 10 5 , 10 6 , 10 7 , and 10 8 , respectively. Sands A and B are added at equal rate tab = 0.5. 



of additions of each sand. Rate equations for the reduced 
variables read 
dX 



dt 

dX A 
dt 

dX B 
dt 

dZ 
dt 

dZ A 
dt 

dZ B 
dt 
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— 2X (Z + Za 
-{&a + Sb)X , 
(X — Xa){Zq + Za 
+5aX - SbXa, 
(X - X B )(Z + Z A + Z B ) + Z B 
+6bX — 5aX b , 

(a A X A + a B X B - 2Z )(Z Q + Z A + Z B ) 
+ol b 5aXb + olaSbXa, 
{(1 - a A )X A + Z }{Z + Z A + Z B ) - Z A 
+ (1 - a A )S B XA, (14) 
{(1 - a B )X B + Z }{Z + Z A + Z B ) - Z B 
+ {l-a B )S A X B (15) 

where aA(t) = X(i t0 \/X A = exp(— S A t) and a B (t) = 
X(oa)/ x b = exp(-5 B t). 

Numerical solutions of the equation system 110|l - l|15|) 
are plotted in Fig. |3J It is found that in the long time 
limit, the system goes to a steady state which is sup- 
posed to be SOC. The steady state is able to be exam- 
ined analytically. When t — > oo, then a.A,ot B — > and 
Xq, Zq — > 0. Therefore the following relations are derived 
from Eqs. (QUJ-CSJl- 

Za Z U (16) 
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FIG. 3: Time evolution of the reduced six variables at 5a 
8b = 0.0001. Initial values are set to Xo = 1.0 and Xa 
X b = Z = Z A = Z b = 0.0. At the steady state, X = Z 
0, X A = X B ^ 0.5, and Z A = Z B ^ V5a/2 = 0.005. 
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we obtain 



Z B > 5a, S b and Z A ■ Z B = 8 A ■ Si 
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(17) 



Za + Z B + 5a 



Equations (|17f) show that the model goes to a critical 
state (Za = Z B = 0) in the limit 6 A , S B ~ 5 — > 0. Out of 
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two conditions for SOC of the first kind, our model 
can successfully remove one of them with respect to the 
control parameter of dissipation e. Furthermore, the or- 
der parameters behave as Za, Zb ~ S 1 ^ 2 when the con- 
trol parameters go to zero 6 — ► 0. These two character- 
istic features of two-component sandpile models, that is 
the existence of the steady state without dissipation and 
the asymptotic behavior to the critical state in 5 — > 0, 
distinctively differ from usual SOC models and are strong 
evidences that the model belongs to different universal- 
ity class. Consequently, we are successfully able to find a 
new mechanism and class of SOC. So, we call this SOC 
of the second kind. 

We have investigated two-component sandpile mod- 
els with the rule (b) and periodic boundary conditions 
mainly on 1— dimensional lattices. Extention to two com- 
ponents is essential to generate an infinite number of 
stable states which substitute for boundary dissipation. 
Therefore, the model becomes more natural because the 
dissipation is usually introduced artificially only for stop- 
ping avahanches. Simulation results indicate the model 
shows SOC behavior. The mean-field theory confirms 
that the model goes to the critical steady state belong- 
ing to a different universality class from one-component 
models. Here, only rough-tuning of rate constants of ad- 
ditions S is enough to induce SOC. We can successfully 
construct a new SOC model without introducing any ar- 
tificial dissipation mechanisms. Consequently, the two 
conditions J5J) are not the necessary condition for SOC, 
which is the answer to the question mentioned previously. 

It is interesting to summarize SOC models from the 
standpoint of the number of conserved quantities. Mod- 
els that have no conserved quantities, such as con- 



tact processes, show power laws strictly on the critical 
point. Therefore, fine-tuning of control parameters is re- 
quired. Models that have one conserved quantity, such as 
one-component sandpile models, become SOC in rough- 
tuning of two control parameters, an addition and a dis- 
sipation rates of sands, around zero. Two-component 
sandpile models with two conserved quantities show SOC 
in rough-tuning of only one control parameter, an addi- 
tion rate of sands, around zero. It could be concluded 
that the more conservation laws result in the less control 
parameters. 

Sandpile models are considered to have some rel- 
evance to power laws of earthequake magnitude dis- 
tributions, called Guthenberg-Richter law[14|. One- 
component sandpile models have to introduce some arti- 
ficial boundary dissipations in order to stop avalanches. 
However, the surface of the earth has no apparent bound- 
ary. At this point, it is more natural to apply two- 
component sandpile models to earthequake magnitude 
distributions. Furthermore, an earthquake could be trig- 
gered by multiple physical quantities, such as elastic 
energy and stress of earth's crust. If an earthquake 
takes place only when both the energy and stress reach 
threshold simultaneously, the distributions of earth- 
quakes would be well described by the SOC of the second 
kind. In this letter, we express the process in terms of a 
sandpile. Obviously, sands are merely symbols and could 
be any kind of substances which trigger the dynamics. In 
addition to earthquakes, therefore, a wide variety of SOC 
phenomena triggered by multiple factors would belong to 
SOC of the second kind. Future studies will make clear 
these points. 
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